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Abstract 

Recently the complete characterization of a general Gaussian dissipative system having a unique 
pure steady state was obtained in [Koga &: Yamamoto 2012, Phys. Rev. A 85, 022103]. This 
result provides a clear guideline for engineering an environment such that the dissipative system 
has a desired pure steady state such as a cluster state. In this paper, we describe the system 
in terms of a quantum stochastic differential equation (QSDE) so that the environment channels 
can be explicitly dealt with. Then a physical meaning of that characterization, which cannot be 
seen without the QSDE representation, is clarified; more specifically, the nullifier dynamics of any 
Gaussian system generating a unique pure steady state is passive. In addition, again based on the 
QSDE framework, we provide a general and practical method to implement a desired dissipative 
Gaussian system, which has a structure of quantum state transfer. 
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I. INTRODUCTION 



Towards quantum state preparation, which clearly plays a key part in quantum infor- 
mation processing, recently several dissipation-based approaches have been proposed. The 
basic idea of those approaches originates from the trivial fact that a thermal environment 
drives any state to the stable ground state. However, it has been shown that we can some- 
times engineer a desired dissipative environment such that the corresponding stable state 
is a nontrivial and useful one, e.g., a highly entangled pure state [lU]. More specifically, 



under some conditions we are allowed to synthesize an open quantum system described by 
the Markovian master equation 
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^ = pj + [l,p,lI - -LlUp, - -p,LlL,) , (1) 
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such that pt must converge into a given desired pure state poo', that is, the Hamiltonian H 
and the dissipative channel Lk {k = 1, . . . ,m) are appropriately synthesized to achieve this 
goal. One of the main advantages of this approach is that the target state poo is clearly 
robust against any perturbation to the state pt during the dynamical process. In particular, 
it is independent on the initial state preparation. 

In the finite dimensional case a necessary and sufficient condition for Eq. ([1]) to have 
a pure steady state was obtained in j^, 0], and especially in 0] the authors provided a 
sufficient condition for poo to be unique. The uniqueness characterization is of particular 
importance, because without such condition the desired convergence into the target state 
cannot be guaranteed. For infinite-dimensional systems, on the other hand, in pLSn] the 
author particularly focused on a general Gaussian dissipative system and provided a complete 
parameterization of the system having a unique pure steady state. The merit of focusing 
on the class of Gaussian systems lies not only in its importance in quantum information 



technologies [15|, Il6[ but also in the fact that the parameterization is obtained in an easily- 
tractable manner in the phase space; actually the uniqueness of poo can be readily checked 
by simply calculating the rank of a specific matrix, while in the finite-dimensional case 0] 
we are required to verify that there is no specific subspace in the Hilbert space. 

In this paper, we study a Gaussian system having a unique pure steady state in terms 



of a quantum stochastic differential equation (QSDE) [17H20|. The use of a QSDE allows 
us to describe the dynamics of an open system in a form where the stochastic environment 
channels appear explicitly. The master equation ([T]) is obtained as a result of averaging out 
all such stochastic effects brought from the environment. 

The contribution of this paper is twofold. The first one is that we clarify the physical 
meaning of the conditions for the Gaussian system to have a unique pure steady state, which 
cannot be clearly seen when dealing with only the master equation. This result is obtained 
through investigating the QSDE of the corresponding (complex) nullifier. In general, it is 
known that any pure Gaussian state can be characterized as the common zero eigenstate 



of the corresponding nullifier operators |21|; this is the reason why investigating the full 
behavior of the nullifier provides new information about the dynamic process towards the 
target pure state. Actually, we show that the nullifier dynamics of any Gaussian system 
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generating a unique pure steady state is passive. As a byproduct, the result is used to show 
a certain trade-off between the closeness of the steady state to a target Gaussian cluster 



22|-|25| and the convergence time into that steady state. 



state 

In the previous result [13[, although the mathematical characterization of the desired 
dissipative channel Lk was obtained, its actual implementation was not discussed. Actually, 
the resulting desired dissipative channels usually have to non-locally act on the system, and 
no general method to effectively implement such dissipative channels is known. The second 
contribution of this paper is to give a partial answer to the question of how to practically 
construct a desired dissipative system. The proposed scheme has a structure of quantum 



state transfer horn light to a matter [26H28J; more specifically, a desired state of light is first 
generated and then that light field interacts with the oscillator system (memory), which as 
a result acquires the desired state by dissipation. This scheme is indeed practical, because, 
as shown in [2^ , any pure Gaussian cluster state of light can be effectively generated using 
some beam splitters and OPOs. Note that the QSDE approach actually has to be taken in 
order to explicitly describe the input light field. 

We use the following notations: for a matrix A = (aij), the symbols A'^ , A'^ , and A" 
represent its Hermitian conjugate, transpose, and elementwise complex conjugate of A, i.e., 
A^ = (a*j), A~^ = (aji), and A^ = (a*^) = {A^)~^ , respectively. For a matrix of operators, 
A = (ciij), we use the same notation, in which case a*j denotes the adjoint to dij. In denotes 
the n X n identity matrix. 3? and 3 denote the real and imaginary parts, respectively. 

II. PRELIMINARIES 

In this section, a brief introduction to a Gaussian system and its QSDE representation is 



given. Then we review the result of 13 



A. Gaussian dissipative systems 

A general n-mode bosonic system consists of n subsystems with canonical conjugate pairs 
(gj, pi). Denote the vector of total system variables by x := (gi, . . . , g„,pi, . . . ypn^ ■ The 
canonical commutation relation [g^, pj] = i6ij then leads to 



xx'^ — (xx'^)'^ = zS, S 



/„ 

-In 



Now let p be the density operator of this system and write the mean vector by (x) and 
the covariance matrix hj V = {AxAx'^ + {AxAx'^)'^) /2, Ax = x — (x), where the mean 
{X) = Tr (Xp) is taken elementwise. Note that the uncertainty relation V + zS/2 > 
holds. A Gaussian state can be characterized by only the mean vector and the covariance 
matrix. A particularly important fact is that the covariance matrix V corresponding to a 



pure Gaussian state always has the following general representation 2l|, |29 



1 -J- / y-1/2 
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where X and Y are n x n real symmetric and real positive definite matrices (i.e., Y = 
Y~^ > 0), respectively. In other words, a pure Gaussian state is completely parameterized 
by X and Y. An important merit of this representation is that the complex graph matrix 
Z := X + iY can be used for a graphical calculus for several Gaussian pure states |21j. In 
particular, a pure Gaussian state \ipz) having the covariance matrix ([2]) always satisfies 



(3) 



where the equation means that each entry of f acts on \ipz)- Conversely, if a pure Gaussian 
state lip) satisfies {—Z, In)x\ip) = 0, then we have lip) = li'z)- The vector of operators r is 
called the nullifier for the pure Gaussian state l^pz)- 

A linear system is such that the Hamiltonian H and k-th dissipative channel in Eq. 
([1]) are respectively characterized by 











r\ipz) = 0, f := {-Z, In)x = 
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where G = G^ G R^"^^" and Ck G C^". For this system, the time-evolution of (xt) and Vt 
with the state pt obeying Eq. ([T]) are given by d{xt)/dt = A{xt) and dVt/dt = AVt+VtA^ +D, 
respectively. Here, A = E[G + Q{C^C)] and D = S3fJ(CtC)S^ with C = (ci, . . .,CmV G 



Qmx n j^ggg j2gj f^j, morc detailed discussion). In this paper, we assume that the initial state 
of the dynamics is Gaussian; then, at any given time t the state is also Gaussian with mean 
(xt) and covariance Vt, hence let us call such a linear system the Gaussian system. A steady 
state of the Gaussian system exists only when A is a Hurwitz matrix, i.e., all the eigenvalues 
of A have negative real parts. If it exists, the mean vector is (xqo) = and the covariance 
matrix V^o is given by the unique solution to the following matrix equation: 



AV^ + V^A^ + D = 0. 



(5) 



B. The QSDE framework 

The situation we have in mind is that the system interacts with countable set of envi- 
ronment channels. The time-evolution of an observable of this open system is described in 
terms of a QSDE. A most simple form of this equation is obtained when the environment 
channels are all independent vacuum fields with ideal Markovian approximation taken. Let 
di{t) be the annihilation operator of the z-th vacuum field; then the Markovian approxi- 
mation means that di{t) instantaneously interacts with the system and satisfies the CCR 
[di{s), d*{t)] = 5ij5{t—s). Define the field annihilation process operator by y4j(t) = Jl^di{s)ds, 
then this CCR leads to the following quantum Ito rule: 

dAidA* = 6ijdt, dAidAj = dA*dA* = dAidA* = 0. (6) 

The system-field coupling in the time interval [t,t + dt) is described by the unitary operation 
U{t + dt,t) = exp[^-{LidA* — L*dAi)], where Lj is the system operator representing the 
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coupling with the z-th vacuum field. Then the system observable at time t, jt{X) = U^XtJt, 
obeys the Ito-type QSDE 

djtiX) = jt{t[H,X] + J2 {L:XL, - -LfUX - -XL*L,)yt 

i=l 

m 

+ {jt{[X,L,])dA*-M[X,L:])dA,), (7) 



i=l 



where an additional system Hamiltonian H has been added. Note that Ut+dt = U{t, t+dt)Uf 
The mean value {jt{X)) is represented using the (unconditional) density operator pt by 
{jt{X)) = Tr {Xpt), which leads to the master equation ([T]). The change of the field operator 
can also be dealt with explicitly; the output field A- := jt{Ai) after the interaction satisfies 

dA', = jtiLi)dt + dAi. (8) 

We are interested in the QSDE whose system Hamiltonian and dissipative channels are 
given by Eq. (jl]). Let us define At = {Ai, . . . , A^)^, then the vector of system quadratures 



xt = (jtiqi), ■ ■ ■ Jt{qn)Jt{Pi), ■ ■ ■ JtiPn)) satisfies the linear QSDE |20|, l30|-|35 



dxt = Axtdt - i^C^At + i^C^ dA\, (9) 

where the system matrices A and C were defined in Section 2 (a). It is easy to see that 
the linear QSDE ([9]) actually leads to the time-evolutions of the mean and the covariance 
matrix: d{xt)/dt = A{xt) and dVt/dt = AVt + VtA~^ + D. Also the output field equation ([8]) 
of A'^ = {A[, . . . , A'^)~^ then becomes 

dA't = Cxtdt + dAt. (10) 



C. The dissipative Gaussian system generating a pure steady state 



In [13|, some conditions for a dissipative Gaussian system to have a unique pure steady 
state were obtained. A particularly useful result from the environment engineering viewpoint 
is the followin g (r ecall Z := X + iY): 



Theorem 1 



131 ]: Let \^ be a given covariance matrix of the form ([2]). Then, this is the 



unique solution of Eq. ^ if and only if the system matrices are represented by 

C = P'^{-Z,h), (11) 
/ XRX + YRY - VY-^X - XY'^V^ -XR + VY'^ \ , . 

-RX + Y-'T^ R J ' ^^^^ 

where P is a complex n x m matrix, i? is a real n x n symmetric matrix, and F is a real 
nxn skew symmetric matrix (i.e., F + F''' = 0), and moreover, P and Q := —iRY — Y~^T"^ 
satisfy the following rank condition: 

rank(P, QP, . . . , Q"~^P) = n. (13) 
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This theorem states that any dissipative hnear system having a unique pure Gaussian 
steady state is completely parameterized by the three matrices P, R, and F, which further 
have to satisfy the rank condition ( !T3|) . In [13|, this result was obtained through a fully alge- 
braic treatment of Eq. ([5]), and the physical meanings of the conditions were not discussed. 
As mentioned in Section 1, nevertheless, they will be clarified within the QSDE framework; 
for convenience of the later discussion, we note that G satisfies 

GS- (-f) = ( ) Q. (14) 



In J \ In 



III. DYNAMICS OF THE NULLIFIER 

As seen in Eq. ([3]), the pure state {ipz) is the common zero-eigenstate of the nuUifier 
vector f = {—Z, In)x. Hence it is worth to see the time-evolution of ft, when the conditions 
shown in Theorem 1 are satisfied. Noting Eqs. (fTTj) and (fT^ . we have 

(-Z, I) A = (-Z, nm + -(-z, i)^(c^c - c^c^) 
= Q^i-zj) + ^(-/, -z){ ( -f ) p«p^(-z,/) - ( -/ ) ppH-zKi)} 

7^-7 

= Q^(-Z, /) + —^pip^{-Z, I) = (gT _ YPip^)i-Z, /), 
{-Z, = -i{-Z, /)S ^ P» = i{Z - Z^)P^ = -2YP\ 



I 

Therefore, multiplying both sides of Eq. ([9]) by (—2', /„) from the left, we have 

dft = (Q^ - YP^P^Ytdt - 2YPUAt. (15) 

Regarding the output process (ITOl) . as Cxt = P^ {—Z, In)xt = P^h, it is written by 

dA[ = P^ndt + dAf (16) 

The coefficient matrix of the dynamics of f has the following property. 

Proposition 2: The matrix — YP^P^ is Hurwitz if and only if the rank condition f[T5]) 
is satisfied. 

Proof: Let b G C" and A G C be the eigenvector and the eigenvalue of — YP^P^ , 
respectively; i.e., {Q^ — YP^P^)b = Xb. Then, multiplying this equation by b'^Y^^ from the 
left, we have 

5t(y-igT _ pttpT)^ _ ^||y-i/2^||2^ 

where || • || denotes the standard Euclidean norm. This immediately yields b^{Q^Y^^ — 
pip^)b = \*\\Y-^/%f. Recall now that Q = -iRY - Y'^T^ , hence QY'^ = 
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—iR — Y ^r^Y ^ is skew Hermitian. Therefore, adding the above two equations yields 
-2\\P^bf = {\ + \*)\\Y-^/%f, and we have K(A) = -\\P^b\\y\\Y~^/%f. Let us here 
assume that P^b = 0. Then, we have Q^b = Xb, and the matrix C := (P, QP, . . . , Q"'~^P) 
satisfies 

f P^b \ j P^b \ 

P^Q^b XP^b 



C'b 



0. 



\ A"-^P ' b ) 



But this is contradiction to the assumption ( IT3l) . thus P 6 7^ 0. As a result, 3fJ(A) is strictly 
negative, implying that the matrix — YP^P^ is Hurwitz. 

On the other hand, if C is not of rank n, there exists an eigenvector of — yP'P"'', say 
60, that satisfies P^6o = 0. Then, from the above discussion, the corresponding eigenvalue 
A satisfies 3fJ(A) = 0, hence - YP^P^ is not Hurwitz. ■ 

Based on the above result, we can verify that the target pure Gaussian state is certainly 
generated. To see this, let us multiply all the entries of the nuUifier dynamics (fT5|) by the 
system-field composite state vector = \ip) ® |0) from the right (|0) is the vacuum state). 
Then, due to the relation dAt\0) = 0, we have 



d_ 

di 



ftl^) = {Q'^ -YP^P^)rt\^). 



The Hurwitz property of the matrix — YP^P^ is equivalent to asymptotic stability of 
the dynamics, thus the nullifier vector f(|\l/) = f/j*ff/(|\I^) converges to zero. Therefore, in the 
Schrodinger picture, Uool"^) is the common zero-eigenstate of f, meaning that the system 
state becomes the target pure Gaussian state {ipz) as t — )■ 00. 

The Hurwitz property obtained above allows us to obtain a specific input-output relation 
from the incoming field At to the outgoing field A't, through Eqs. (fT5l) and (fT6l) . For 
this purpose, it is convenient to move into the frequency domain where both of these field 
operators as well as the internal system variable Xt are all Fourier transformed. The Fourier 
transformation of Eqs. (fTSi) and ( fT6l) are given by 



iuf{u) = (g^ - YP^P^ 
A'iu) = P'^fiuj)+Aiu) 



2YP^A{ 



00) 



where the tilde notation denotes the Fourier transformed operator. Note that for instance 
Ai{uj) is not the Fourier transformation of but di{t). Also, more precisely, we should 
take Laplace transformation aj(t) — )■ Ai{s) and set s = +0 + iu to obtain the Fourier 

As a result, we have the input- 



transformation; for the rigorous treatment, see |3Ch32 
output map from A{u) to A'{uj): 



A'{uj) = F{uj)A{co), F{co) := - 2P'^ {icu - + YP^P'^y^YPK 



(17) 



Here we have used the Hurwitz property to justify that the initial contribution of the system 
was ignored in Eq. f|T7|) . The m x m matrix F{u)), called the transfer function matrix, has 
a striking property as shown below. 
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Proposition 3: The transfer function matrix F{uj) is unitary for all u. 

Proof. The proof directly follows from Lemma 2 of j3o| , but here it is given for convenience 
of readers. First, to simplify the calculation, let us define P = Y^^^P and Q = Y^I'^QY~^I'^ . 
As Q = —iRY — Y~^T~^ , Q is skew Hermitian; Q + = 0. With this notation, the transfer 
function matrix is represented by F^u) := Im — 2P~^ {iu — + P^P^)~^PK Therefore, we 
have 

F{uj)^F{uj) = im- 2P^{-tu - g» + pip^y^pi - 2P^{iu - + p^p^y^p^ 

+ 2P^{-iu -Q^ + P^P'^)~\2P^P'^){iu -Q^ + P^P^Y^PK (18) 
Here it follows from Q + Q"^ = Q that 

2pflpT ^ - gtt + P«P^) + {iuo -qY + P^P^). 

Substituting this expression for the last term of Eq. (ITSl) . we end up with the relation 
F{uyF{u) = Im, hence F{co) is unitary for all u. ■ 

This result states that the output power spectrum is flat in all the frequency domain, 
i.e., {A'{u)A'{uy) = F{u){A{uj)A{ujy)F{ujy = /„. This means that, at steady state, it 
is impossible to extract any information about the internal system, as long as the matrix 
g^ — YP^P^ is Hurwitz and the input fields are in vacuum or coherent states. 

Now we arrive at the stage where the physical meanings of the conditions given in Theo- 
rem 1 can be clarified. First of all, the structure of Eqs. f|T5|) and f|T6|) as well as the above 
two propositions remind us that the dynamics of the nuUifier is a generalization of that for 
the simple single-mode optical damped cavity whose QSDE is described by 

dat = [iA — —^dtdt — ^JTidAt., dA[ = y/ndtdt + dAt, (19) 

where dt and At denote the intra-cavity mode and the incoming vacuum field mode, respec- 
tively. A and /t denote the detuning and the damping rate, respectively. Clearly, the state 
evolves into the vacuum, and also we have {A' {uo)* A' {u)) = {A{uj)* A{uj)) = 1 for all u, as 
in the nullifier case. These properties arise due to (i) that energy is not supplied through 
the Hamiltonian, (ii) that the field does not supply energy but simply brings about the 
damping of the system, and (iii) that the system is asymptotically stable. Mathematically, 
the first two statements mean that the dynamics does not contain the creation operators 
d^ and A^. Actually, regarding the first one, if the cavity contains a degenerate parametric 
amplifier, which is described by the Hamiltonian Hdpa = i{d*'^ — a^), then the QSDE needs 
to be described in terms of both d and d*. The last condition (iii) guarantees that the state 
uniquely converges into the vacuum as well as that the output field does not contain any 
information about the system at steady state. Systems having the properties (i)-(iii) are in 



general called passive systems [30|, I33M35 



The above discussion implies that the nullifier dynamics is passive; more precisely, we 
obtain the physical meanings of the conditions (ITTl) . (|T2l) . and (fT3l) as follows. 



The matrix C has the form given in Eq. ffTTj) so that the creation process Ai does not 
appear in the QSDE of f^; as mentioned above, this is equivalent to that there is no 
energy supply from the environment to the nullifier. 
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The matrix G has the form given in Eq. (fT2|) so that the corresponding Hamiltonian 
H = x~^Gx/2 does not supply energy for the nuUifier. 

The rank condition ( fT3i) implies that the coefficient matrix — YP^P"^ is Hurwitz, or 
equivalently the asymptotic stability of the dynamics of the nullifier. This guarantees 
that the output power spectrum is fiat in all frequencies, meaning that the output 
field does not contain any information about the system at steady state. 



The last statement can be understood by studying the filtering equation [18|, |36[, which 
enables us to update the conditional state based on the measurement result of the output 
field. In general, when the state of the master equation reaches the steady state and it is 
pure, then the corresponding filtering equation is identical to the master equation, meaning 
that we do not obtain any new information through measuring the output field for updating 
our knowledge. Note that this does not mean that the system is not controllable. 



IV. QUANTUM STATE TRANSFER FOR THE DISSIPATIVE SYSTEM ENGI- 
NEERING 

We have seen in Theorem 1 how the dissipative channel Lk = cjx with C = {cj , . . . , c^) 
should be chosen to engineer a desired Gaussian dissipative system. When we aim to gener- 
ate a certain (useful) Gaussian state, however, it often turns out that the resulting Lk has to 
possess a specific structure which is hard to actually implement. For instance, a dissipative 
channel interacting with all the nodes, i.e., L = ^igi +£2'?2 + • • • + C^n, will be hard to imple- 
ment. In this section, for the specific case where the system is subjected to n independent 
input optical fields (i.e., m = n), we provide a practical procedure for implementing desired 
dissipative channels. 

Let us first introduce the field quadratures {k = 1, . . . ,n) 

Qk = iAk + Al)/V2, Pk = {Ak-Al)/V2i, (20) 

which satisfy the canonical commutation relation [dQi, dPj] = 6ijdt. Then, defining Q = 
(Qi, . . . , QnV and V = {Pi, . . . , PnY , we find that the QSDE @i is rewritten by 



dxt = Axtdt + y2S(C7, Cj)dWt, = , (21) 




where = ^{C) and d = Q{C). 

The situation we have in mind is that a desired pure Gaussian state of light is first gen- 
erated, and then, that state is transferred to the system through the system-field coupling; 



see Figure 1. This is the framework of the quantum state transfer |26l-l28| : in this case the 
system is called the memory and it should be independent on the input state we will trans- 
fer. More specifically, the target mode Wf is obtained from the vacuum mode Wt through 
the transformation = SVVt where the symplectic matrix 5* is given in Eq. (j2]); then the 
quantum Ito rule gives 

dVV^idW^y = SdWtdWjS^ = S{In + iT.)S'^dt/2 = {SS^/2 + zS/2)cit, 
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FIG. 1: Quantum state transfer from the input mode to the memory mode xt. 



implying that the covariance matrix of the input field VVf is certainly V = SS~^ /2 (in the 
rigorous sense this statement should be given in terms of the power spectrum density; see 
j31|). Note that 5* can contain a squeezing process, which was not included in the original 



QSDE framework of Hudson and Parthasarathy [l7|; see [32] for a detailed discussion. Now 
the system (!2T|) is written as 

dxt = Axtdt + BdWf, B := V2^{Cj , Cj)S-\ (22) 

where, as shown above, the new input field VVf carries information of the target Gaussian 
state. The system, which serves as a memory, should satisfy the following two requirements: 

(Rl) The memory system ( 122|) should not possess any information about the input state; 
that is, the system's coefficient matrices A and B should be independent on Z = 
X + iY. 

(R2) The state of the memory system ( |22|) should converge to the target Gaussian state 
with covariance matrix Voo = SS~^ /2. That is, C and G should be of the form ( fTTjl 
and ( fT2l) with P and Q satisfying the rank condition ( fT3l) . 

Below we give a characterization of the desired memory system: 

Proposition 4'- Assume that the system satisfies the requirements (Rl) and (R2). Then, 
the system has to be of the form 

dxt = -2K^xtdt - 2Krf>Vf , (23) 

where k is a scalar constant. 

Proof. Note that ^"^ = S^^S^. Then, substituting C = P^{-Z,In) for B in Eq. ([22]), 
we have 

-P2Y'/' - P^Y-^/'X P,Y-^I^ 
^ V -(^^1 + XP2)Y^I^ - (XPi - YP2)Y-^/^X {XPi - YP2)Y~^I^ 

where Pi = 3fJ(P) and P2 = Q'(P). First let us look at the (2,2) block matrix; since X can 
take any symmetric matrix, here we set X = 0, implying YP2Y~^^'^ is independent on Y . 
This readily implies that P2 has to be of the form P2 = \/2kY'^^/'^ with k a constant. Then, 
XPiF~^/^ has to be independent on X and Y . But as the (1,2) block matrix 9 := PiF~^/^ 
also has to be independent on X and Y , thus this is the case for XQ as well. Then, B = 



10 



is only allowed, hence we obtain Pi = 0. With these selection of Pi and P2, the (1,1) and 
(2,1) block matrices of B take — 2k/„ and zero, respectively. As a result, B = —2nl2n- 

Next let us consider the matrix A = S(G + '^{C^C)). From the above discussion, now we 
have C = \/2kY~^^'^{—Z, In), which leads to A = T^G — 2/t^/2n- This means that the matrix 
G given in Eq. f[T^ must be independent on X and Y. Then, similar to the above discussion, 
by setting X = 0, we find that G = {YRY, TY'^ ; F-^F^, R) has to be independent on X 
and Y. But this requirement is only satisfied when R = and F = 0; as a result, we have 
G = 0. m 

This proposition states that, in order to dissipatively generate a desired pure Gaussian 
state in the state transfer setup, we are required to prepare identical and independent 
oscillators as memories. Note that any pure Gaussian cluster state (see the next section) 
can be effectively generated from the vacuum fields by applying suitably combined two-mode 



squeezing Hamiltonians and beam splitters 25|, hence the proposed scheme is practical. 



V. EXAMPLES 



A. Example 1: Gaussian cluster state generation 



It was shown in [21|] that the graph matrix Z = X + iY can be used to capture sev- 
eral Gaussian graph states in a convenient graphical manner. In particular, the so-called 
canonical Gaussian cluster state 22l-l25| , which plays an essential role in continuous- variable 
one-way quantum computation, corresponds to 



Z = X + ie-"^"!, 



(24) 



where X is the symmetric adjacency matrix representing the graph structure of the cluster 
state: for instance, the matrices 



X 



/O 1 0\ 
10 10 
10 1 

Vo 1 0/ 



X 



/O 1 1 1\ 

10 

10 

V 1 / 



X 



/O 1 1\ 
10 10 
10 1 

V 1 1 0/ 



represent the chain, T-shape, and square structures, respectively (see Figure 2). On the 
other hand, Y = e~^°J corresponds to the approximation error of the state with covariance 
matrix ([2]) to the ideal cluster state; that is, bigger a means that the state well approximates 
the ideal cluster state having the graph structure assigned by X. 
Now, the nuUifier dynamics f[T^ is of the form 



dft = (g^ - e~^''P^P^)rtdt - 2e~2ap«^^^ 



(25) 



and the real part of the eigenvalue of the coefficient matrix — e ^apflpT s)fj(^^^ _ 
P^6||, where b is the corresponding eigenvector. This means that making a bigger. 



— e" 



-2a I 



or equivalently making the state more close to the ideal cluster state, renders the stability 
of the nuUifier dynamics worse. Another observation from a more practical viewpoint is as 
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FIG. 2: Typical cluster state with (a) chain, (b) T-shape, and (c) square structures. 

follows; let us define the convergence time to the target by T = l/min|3?(A)| and denote 
the approximation error of the state to the ideal cluster one by e = e~^". Then, it is 
straightforward to find Te > c with c a constant. Therefore, in order to dissipatively 
generate a pure Gaussian state that is very close to a desired cluster state, the convergence 
time has to be long. 



More generally, as shown in [21|, the matrices X and Y respectively correspond to the 
ideal and realistic parts of a Gaussian graph state in the sense that the covariance matrix 
of {—X,In)x is given by Y/2. That is, Y can be regarded as the approximation error in 
approximating the ideal graph structure X. Therefore, the above-mentioned trade-off holds 
for a general Gaussian graph state. ■ 



B. Example 2: Two-mode squeezed state 

There exist a number of proposals to generate a steady two-mode squeezed state in for 
instance atomic ensembles or nano-mechanical oscillators. The system matrices describing 
the two-mode squeezed state are given by X = and 



Y 



cosh(2a) — sinh(2a) 
— sinh(2a) cosh(2a) 



where a denotes the squeezing parameter representing the degree of entanglement. In [l^ 
the dissipative channels achieving this goal was shown to be 

Li = /ifli + z/flg, L2 = fia2 + i^al, (26) 

where n = cosh(Q;) and u = — sinh(a;), while H = 0. In our formulation, this corresponds 
to setting (lit : 

/ icosh(a) isinh(a)\ 



isinh(a) icosh(a) 

However, the dissipative channels ( l26l) are not easy to implement, since they are global and 
nontrivial coupling between the systems and the environment. 

On the other hand, this dissipative system can be more easily implemented within the 
state transfer framework provided in Section 4, because we are only required to generate a 
two-mode squeezed state of optical fields and prepare two identical and independent oscil- 
lators. Note that a two-mode squeezed state of light can be effectively generated using a 
non-degenerate OPO. ■ 
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VI. CONCLUSION 



In this paper, the dissipation-based state preparation method for general Gaussian case, 



which was originally formulated in [13| , was reconsidered in terms of the QSDE. This ap- 
proach clarified that the nullifier dynamics of any Gaussian system generating a unique pure 
steady state is passive. As a byproduct, it was shown that there exists a trade-off between 
the closeness of the steady state to a given ideal graph state and the convergence time to 
that state. In addition, a convenient physical implementation method of a desired Gaussian 
dissipative system was provided; the scheme has the structure of quantum state transfer, 
which is a key ingredient in quantum information technologies. 
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